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Abstract 

We introduce a magnetic hard hexagon model with two-body restrictions 
for configurations of hard hexagons and investigate its critical behavior by 
using Monte Carlo simulations and a finite size scaling method for discreate 
values of activity. It turns out that the restrictions bring about a critical 
phase which the usual hard hexagon model does not have. An upper and a 
lower critical value of the discrete activity for the critical phase of the newly 
proposed model are estimated as 4 and 6, respectively. 

PACS number: 64.60.-i, 75.10.-b 

I. INTRODUCTION 

Critical behavior of the hard hexagon model on a triangular lattice is well known through 
its exact solution by using the corner transfer matrix@. Below a critical value of the activity 
z ~ 11, the system is in the disordered phase. Therefore, a correlation function between 
hard hexagons decays exponentially there!. Above the critical value of the activity, the 
hexagons on the lattice condense and make an ordered configuration. The degeneracy of the 
ordered configuration is triplefold and the universality class of the critical point is the same 
as that of the three-state Potts modelBH. The fact that the hexagons are hard means that 
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there is a restriction between the nearest-neighbor sites, such that, any two hexagons can 
not occupy neighboring two sites simultaneously on the triangular lattice. This property of 
the hard hexagon model can be regarded as the nearest-neighbor exclusion. Adding to the 
nearest-neighbor exclusion, a next-nearest-neighbor exclusion causes fourfold degeneracy in 
the ground state. Todo and SuzukiB claimed that the hard hexagon model with the next- 
nearest-neighbor exclusion belongs to the same universality class as that of the four-state 
Potts model by means of the numerical study by using the phenomenological renormalization 
group method and the conformal invariance. 

Recently, it was found that an antiferromagnetic spin-,? Ising model on the triangular 
lattice has a relation with the hard hexagon model with many restrictions for configurations 
of hexagons!. The phase transition as a function of S in the ground state of the spin-S" Ising 
model was studied by using a mapping from the ground state degeneracy of the system at 
zero temperature to a partition function of a spin-1/2 system at a pseudo temperature!!; 
the mapped system is called a A model. The original 5=1/2 Ising model was exactly 
solved and it is well known that the system is critical with 77 = 1/2 at zero temperature!!. 
Nagai et al.i showed that the critical index 77 decreases from 1/2 for S = 1/2 when the 
spin S increases from S — 1/2. It has been clarified that 77 becomes zero at S = 7/2 
discontinuouslyil. This behavior of 77 suggests an appearance of a partial long-rang order for 
the system with S > 7/2; it was shown explicitly for S = oo0. We remind that the A model 
is the hard hexagon model with many kinds of restrictions!. In contrast with the usual hard 
hexagon model which does not have any extra restrictions but only the nearest-neighbor 
repulsion, the A model is critical for 1/2 < S < 3!. 

The purpose of the present paper is to introduce a new hard hexagon model; we call 
it a magnetic hard hexagon model. The magnetic hard hexagon model with up to the 
third-nearest-neighbor repulsions is investigated by Monte Carlo simulations and a finite 
size scaling for evaluation of 77. We find that the new hard hexagon model has an ordered 
phase, a critical phase and the disordered phase. We discuss a relation of the new hard 
hexagon model with the spin-S" Ising model on the antiferromagnetic triangular lattice. 



In section ||, we introduce the magnetic hard hexagon model with exclusions up to the 
third-nearest-neighbor sites with relation to the spin-5 1 Ising model on the triangular lattice. 
In section [TTI|, we show our results about a phase diagram of our new model and discuss 
these results. Concluding remarks are given in section |TV]. 

II. MAGNETIC HARD HEXAGON MODEL WITH RESTRICTIONS 

Let us introduce a magnetic hard hexagon model which is related to a representation of 
the ground state degeneracy for the spin-S 1 Ising model on an antiferromagnetic triangular 
lattice (ATL). The "occupation number" of magnetic hard hexagons is denoted by Q at a 
site i on the triangular lattice A; Q takes on {— 1,0, +1}. If the site i is occupied by an 
" up magnetic hard hexagon" or by a " down magnetic hard hexagon" , then we assign them 
d = +1 and —1, respectively. If the site i is empty, then we set Q = 0. 

The ground partition function of the system is written as 

^RMHH = E^ l6l II( 1 -|Ci|lOl) 

{6} <M> 

x n a*) n ki-d(i) (i) 

(i,ky z (i,iy A 

where N = |A|, and the products with (i,k)' and are taken over the nearest- 

neighbor pairs of sites, the next-nearest-neighbor pairs of sites and third-nearest-neighbor 
pairs of sites, respectively. We explain below a relation of Zrmhh with the ground state 
degeneracy of the spin-.S' Ising model on the ATL. In this relation, the activity denoted by 
z is equal to 2S. 

The Hamiltonian of the spin-S* Ising model on the ATL is written as: 

K{S} = J y £S i S j , (2) 

where J(> 0) is the interaction constant between nearest-neighbor sites denoted by 
and Si is a spin variable which takes a value on {—S, —S + 1, • • • , S — 1, S}. It has been 
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shown that the ground state degeneracy of this system is equivalent to a partition function 
of a spin- 1/2 Ising model as follows!!: 

Z= ]T exp{ln(2S)^{a,}} , (3) 

{<7 i |i?{ CTi }/L2=-J/4} 

where the summation for spin configurations {07} is taken over the ground-states configura- 
tions of spin-1/2 Ising model on the ATL with 07 G { — 1/2, +1/2}: the energy of the system, 
E{ai}/L 2 , has to be — J/4. We recall that this is the partition function of the A mo del. 
The linear size of the system is denoted by L. K{a{\ denotes the total number of free spins 
and is expressed as follows: 



K{a l } = Y,k i {ai} , (4) 

ieA 




where in Eq.(|]) the summation is over the set of sites in the lattice A and in Eq.fl5|) C\(i) 
and C 2 (i) are sets of sites, each of them composed of three sites encircling the site i = (x, y) 
as shown in Fig.|I[ When right bonds shown by thick lines are arranged in such a way given 
in Fig||, the value of the function ^{07} becomes 1; this means that a hexagon occupies the 
site i = (x,y). In this situation, the site (x + l,y) cannot be occupied by another hexagon, 
because the value of V(x,y)V{x+\,y) on the bond (x, y) — (x + 1, y) is —1/4 and hence the value 
of k( x+ i^){ai} is 0. One cannot put another hexagon on the site (x, y — 1) because the value 
£7 (x,»)°'(a!,j/-i) is +1/4. By the same reasons as these for the sites (x + l,y) and (x,y — 1) 
explained above, all of sites encircling the site i = (x, y) cannot be occupied by another 
hexagon when the site % is occupied by a hexagon. Therefore, we can regard ki{a{\ as an 
occupation number of the hard hexagon. 

In the partition function written in Eq.(^), the summation is taken over the ground state 
spin configurations. If we regard the function ki{a{\ of the spin configuration {07} as a 
variable hi, the partition function is rewritten as a summation over configurations of hard 
hexagons {ki} as follows: 
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Z = $>xp{ln(2S)$>] 

{ki} I ieA ) 



x J2 S(E{a t }/L 2 + J/4) J] S(ki - k^}) . (6) 

{a,} ieA 

where we define a function 5(x) as 

| 1 (x = 0) 

= { (7) 
[ (a; / 0) . 

The restriction for occupations of the nearest-neighbor sites, namely the nearest-neighbor 
exclusion, is included in the expression llieA'K^i — h{ai}) in Eq.(H). One can write it 
explicitly as follows: 



x J2 5(E{a t }/L 2 + J/4) J] 5(ki - hWt}) . (8) 
{^} «eA 

Comparing this form with the usual hard hexagon model with only the nearest-neighbor 
exclusion0, one notices that 25 plays a role of activity. When the partition function is 
expressed in terms of the occupation number of the hard hexagons, we can see that there 
are many other restrictions for configurations of hard hexagons than that for the nearest- 
neighbor sites. In order to embody some of these restrictions, we introduce a magnetic hard 
hexagon in the followings. 

Let us define the "occupation number" of the magnetic hard hexagon Ci{ a i} as follows: 

( i {a l } = 2a i n (^ + 2<W) II (I ~ 2 Wj) 

jec\(t) vz J jec 2 (i) vz J 

= 2a t k { Wi} ■ (9) 

By introducing the function d{(Ti}, one can distinguish two types of hard hexagons; one of 
them is given by cr,j = +1/2 and the other by <7j = —1/2. They correspond to Q = +1 and 
Q = —1, respectively; we call them an up magnetic hard hexagon for Q — +1 and a down 
magnetic hard hexagon for Q = —1. Regarding the occupation number of the magnetic hard 
hexagon as a variable similar to that in the usual hard hexagon, the partition function (^|) 
can be rewritten by a sum over configurations of magnetic hard hexagons {Q} as follows: 




{Cd t > (hi) 



x £ 5(£{a z }/L 2 + J/4) J] *(Ci - d{ai}) . 



(10) 




In this description, the factor ILeA <K0 — Ol ";}) gives extra exclusion conditions as for con- 
figurations of magnetic hard hexagons other than the nearest-neighbor exclusion. Namely, 
the magnetic hard hexagons have in general next-nearest-neighbor exclusion, third-nearest- 
neighbor exclusion and so on, in addition to the nearest-neighbor exclusion. 

Let us assume that a site i = (x, y) is occupied by an up magnetic hard hexagon with 
C(x,y) = +1 as shown in Fig.^. Since two sites (x + l,y) and (x, y — 1) at the boundary 
between two hexagons located at (x, y) and (x+ 1, y — 1) are connected by a right bond with 
&(x+i,y) = — 1/2 and a( Xj y-i) = +1/2, an up magnetic hard hexagon with C(x+i, y -i) = +1 
on the site at {x + l,y — 1) can exist with that with C(x,y) — 1- However, following to the 
definition of magnetic hard hexagon defined by Eq.(^), the site at {x + l,y — 1) cannot 
be occupied by a down magnetic hard hexagon with (( x +i, y -i) = — 1- The other five next- 
nearest-neighbor sites of (x, y) are in a similar situation. Thus we realize that there is a 
restriction between two magnetic hard hexagons on the next-nearest-neighbor sites. 

We notice that there is another restriction due to the definition of magnetic hard hexagon 
as explained in Fig||. Assume that a site at (x, y) is occupied by an up magnetic hard 
hexagon. When we put a magnetic hard hexagon at the site (x + 2, y) for example, this 
magnetic hard hexagon is in contact with the magnetic hard hexagon at (x, y) on the site 
(x + l,y). Since C(x,y) — +1> the value of spin at (x + l,y) is —1/2. When the magnetic 
hard hexagon at (x + 2, y) has C{x+2,y) — 1> the bond between sites (x + 1, y) and (x + 2, y) is 
not satisfied and hence the magnetic hard hexagons at the site (x, y) and (x + 2, y) are not 
compatible. On the other hand, the magnetic hard hexagon with G x +2,y) = — 1 is compatible 
with that with Gx,y) = +1- Other third- nearest-neighbor sites are in the same situation as 
that of the site (x + 2, y). Thus we have the restriction that a magnetic hard hexagon located 
at the third-nearest-neighbor sites, /, of an occupied site i, Q, has to have the reverse sign 
of that of Ci5 = — sign(Ci)- 
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We have explained three kinds of restrictions for configurations of magnetic hard 
hexagons. Although they originally stem from a definition of the magnetic hard hexagon as 
a function of spin configurations defined in Eq.(|9|), we can write them explicitly in terms of 
{Q} as follows: 

R{&} = U0- - MQ) 

x n l^+dCk) n l^-dCi) en) 

(i,k)' 

where (i, k)' and (i, I)" mean the next-nearest-neighbor and the third-nearest-neighbor pairs, 
respectively. In the A model for which its partition function is given by Eq.(TjI]), there are 
many other restrictions than those two-body exclusions mentioned above. It is expected 
that those two-body exclusions give an important effect for a large value of the activity and 
hence the magnetic hard hexagons condense and make an ordered phase for large values of 
activity; the large value of the activity corresponds to a large value of the spin. 

We now propose a new magnetic hard hexagon model which has only three kinds of the 
two-body exclusions and investigate a critical behavior of the magnetic hard hexagon model; 
the partition function of the new model is described as follows: 

z RMH H = E ex p{ in (2S)$: iCiiWo} , (12) 

{&} I «GA J 

where R{(i} is the restrictions defined in Eq.flTTD; this is equivalent to that given by Eq.(^). 
Although the usual hard hexagon model with only the nearest-neighbor exclusion was solved 
exactly, it seems difficult to obtain an exact solution for the present new model. We therefore 
carry out Monte Carlo (MC) simulations for the model. 

III. RESULTS FROM MONTE CARLO SIMULATION AND DISCUSSIONS 

We perform MC simulations at a finite temperature given by T = 1/ ln(2S) with a unit of 
the Boltzmann constant &b = 1- We used the Metropolis dynamics to update configurations 
of the magnetic hard hexagon in the MC simulations. A finite size scaling methodlHll is 
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used to estimate a critical index r] which describes a decay of correlation functions between 
hexagons separated by a long distance. The index 77 is determined from 



where A, B and C mean three kinds of sublattice of the triangular lattice. We need two 
finite lattices at least to obtain the index 77 from this relation as follows: 



In the present study, we set Li, L 2 , ■ ■ ■ , Lq equal to 24, 36, 48, 60, 90, 120, respectively, in 
this order. This method is based on the assumption that a system is in a critical phase. If 
the system is in the disordered phase, 77 estimated by this method would appear to be 2. 
This value does not mean that the system is in a critical phase with the index 77 = 2 but in 
the disordered phase. On the other hand, if the system is in an ordered phase, it should be 
zero. 

In the present MC simulations, we used more than 10 5 MC step per site for each system 
size. Figure |5] shows results of size dependences of critical index 77 for z = 1,2,3 and 3.5. 
We observed that r\ approaches to 2 in the thermodynamic limit for these values of activity. 
These results mean that the systems with these values of activity are in the disordered phase. 
For a small value of the activity, the density of magnetic hard hexagons is low and hence 
there are many lattice sites unoccupied by a magnetic hard hexagon. 

In Figs.^| and [F], size dependences of the critical index 77 are shown for z = 4, 5, 5.5, 6, 6.5 
and 7. For z c \ < z < z c2 , the values of 77 take finite values in the thermodynamic limit 
where 3.5 < z& < 4 and 6 < z c2 < 6.5. Since the values of 77 are obviously less than 2 
and actually less than 1/2, the systems with these values of activity are in a critical phase. 
In contrast to the case with z < z c \, two-body restrictions considered in our model play a 
role of providing a critical phase for z& < z < z c2 . For z = 6.5 and 7, we found that the 
value of 77 becomes zero in the thermodynamic limit. These results suggest an appearance 




(13) 



v{L i+ i) = 2 



ln{A 2 (L m )M z (L,)} 
ln(L i+ i/Lj) 



(14) 
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of ordered phase at z = z C 2- From the relation of the present model with the spin-S" Ising 
model on the ATL, we have z = 2S and hence S c i = 2 and S C 2 = 3. Note that the discrete 
value of activity, namely the value of spin where r] becomes zero is very close to that in the 
original spin system; S C 2 = <M. This agreement, we think, is due to a dominant role of the 
two-body restrictions in the region of large activity, which gives high density of magnetic 
hard hexagons. 

In Fig.|| we show a comparison of phase diagrams for the usual hard hexagon model, 
the magnetic hard hexagon model with two-body restrictions and the spin-S" Ising model on 
the ATL. The usual hard hexagon model with only the nearest-neighbor exclusion has two 
phase, that is, the disordered phase and the ordered phase. There is no critical phase in the 
usual hard hexagon model. On the other hand, the spin- S Ising model on the ATL, which is 
equivalent to a magnetic hard hexagon model with many restrictions, has the critical phase 
and the ordered phase. There is no disordered phase in the original spin-S" Ising model. Our 
new model investigated in the present paper has three phases, that is, the disordered phase, 
the critical phase and the ordered phase. The ground state degeneracy is sixfold. Hence we 
expect that our new model belongs to the same universality class as that of six clock model. 

IV. CONCLUSIONS 

We have proposed a new model: a magnetic hard hexagon model with two-body ex- 
clusions. We have shown a relation of the new model and the spin-5 1 Ising model on the 
antiferromagnetic triangular lattice. The phase diagram of the new model has investigated 
by evaluating the critical exponent 77 by means of Monte Carlo simulations. It turned out 
that the model has three phases, that is to say, the disordered phase for z < z c i, the crit- 
ical phase for z& < z < z c2 and the ordered phase for z > z C 2, where 3.5 < z c \ < 4 and 
6 < z C 2 < 6.5. We notice that an upper critical value S C 2 for the critical phase in the model 
is the same as the critical value of spin in the original spin-5 1 Ising model. 

Although the two-body restrictions considered in our new model are not sufficiently many 
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restrictions when they are compared to those included in the original spin-S* system, namely 
the A model, we saw they are enough to provide a critical phase, which does not exist in the 
usual hard hexagon model. We estimated in the present paper only the index 77 to clarify the 
phase diagram of our new model. Calculations of the other critical indices and estimations 
of more definite values of z ci and z c2 are left as a future problem to understand the critical 
property of our new model. 
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FIG. 1. Set of sites, C\{%) and Ciii), encircling the site i. These appear in Eq.([|) for evaluation 
of the occupation number ki{ai} of hexagon. 




FIG. 2. Restriction between two hexagons on the nearest-neighbor sites. Thick solid lines 
encircling the site %= (x,y) mean right bonds which give ki{a{\ = 1 corresponding to an occupation 
of this site by a hexagon. One cannot put another hexagon, for example, on (x + 1, y) or (x, y — 1) 
obviously. 
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FIG. 3. Restriction between two magnetic hard hexagons on the next-nearest-neighbor sites. A 
magnetic hard hexagon Q = +1 occupies the site i = (x, y), and therefore a magnetic hard hexagon 
(it = +1 can occupy but Q = —1 cannot the site i' = (x + 1, y — 1). 
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FIG. 4. Restriction between two magnetic hexagons on third-nearest-neighbor sites. 
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FIG. 5. Size dependences of critical index rj for z = 1, 2 and 3. The values of rj become 2 in the 
thermodynamic limit. This means that the system is in the disordered phase. 
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FIG. 6. Size dependences of critical index rj for z = 4, 5 and 5.5. The values of rj apparently 
become finite less than 1/2 in the thermodynamic limit. These results mean the systems with 
z = 4, 5 and 5.5 are in a critical phase. 
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FIG. 7. Size dependences of critical index r\ for z = 6, 6.5 and 7. For z = 6, the value of 
r\ becomes finite in the thermodynamic limit. On the other hand, for z = 6.5, 7 the value of r/ 
becomes zero, which corresponds a ordered phase. 
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FIG. 8. (a) Phase diagram of the usual hard hexagon model with only the nearest-neighbor 
exclusion. There is no critical phase, (b) Phase diagram of our model investigated in the present 
paper. There appear a critical phase between the disordered phase and the ordered phase, (c) 
Phase diagram of the spin-5 Ising model on the ATL. There is no disordered phase. 
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